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INTRODUCTION 


Ultrasonic  velocities  have  been  used  to  study  residual  stresses  in  metals 
for  many  years.  The  linear  relationship  betMeen  these  two  quantities,  the 
stress  acoustic  constant,  can  be  expressed  in  terms  of  second-  and  third-order 
elastic  constants.  The  problem  is  that  the  sound  velocity  of  a  metal  is  also 
affected  by  plastic  deformation  and  crystalline  orientation.  If  one  deals  with 
a  polycrystalline  material  of  sufficiently  small  grain  size  with  random  orien¬ 
tation,  the  velocity  measurement  averages  the  dependence  on  orientation.  Thus, 
one  is  still  left  with  the  effects  of  plastic  deformation. 

Plastic  deformation  can  change  the  sound  velocity  and  the  stress  acoustic 
constant.  In  a  practical  situation,  when  one  tries  to  measure  residual  stresses 
via  shear  waves,  the  effect  of  plastic  deformation  cannot  be  taken  into  account 
because  the  amount  of  plastic  deformation  is  not  known  a  priori,  and  also 
because  the  change  of  velocity  with  plastic  deformation  is  not  available.  In 
addition,  one  does  not  know  how  these  quantities  are  affected  by  plastic  defor¬ 
mation  in  compression  versus  plastic  deformation  in  tension. 

In  this  work,  we  have  studied  the  acoustoelastic  constants  for  shear  and 
longitudinal  waves  obtained  in  tension  and  compression,  and  the  effect  of 
plastic  deformation  in  tension  and  compression  on  these  constants.  We  have  also 
investigated  the  effect  of  plastic  deformation  on  the  velocities.  The  specimen 
used  was  a  bar  of  ASTM  4340  steel  which  had  been  hot-rolled.  The  bar  was 
austenitized,  quenched,  and  tempered  to  a  hardness  of  Rc  31.  The  measurements 
were  accomplished  by  continuous  monitoring  of  ultrasonic  time  delays  and  of 
transverse  and  longitudinal  strains  during  the  four-point  bending  experiment 
through  elastic  and  plastic  regions  of  deformation. 
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THEORETICAL  CONSIDERATIONS 

Elastic  ind  Plastic  Deforwation  of  Rectangular  Bars 

Assuming  uniaxial  stress  oq  and  ideal  elastic-plastic  behavior,  the  yield 
condition  for  a  narrow  curved  beam  is  <7g  »  ±  Oq,  with  Cq  taken  as  the  yield 
strength  in  compression  and  tension.  Yielding  begins  at  the  extreme  fibers  y  = 

±  h/2,  where  y  is  the  distance  from  the  neutral  axis,  and  h  is  the  height  of  the 
beam.  With  increasing  load,  the  elastic-plastic  interface  moves  inward.  The 
central  portion  of  the  beam  -p  <  y  <  n  remains  in  the  elastic  state  with  the 
stress  given  by 

<^0  •  Oq  y/n  *  Ey/Ro  (D 

where  E  is  the  modulus  of  elasticity  and  Rq  is  the  radius  of  curvature  of  the 
beam.  Stresses  in  the  two  plastic  regions  are  given  by  s  ±  Oq  and  independ¬ 
ent  of  the  strain  €  ■  y/Ro. 

During  unloading,  the  plastic  region  reverts  to  the  elastic  state  and  the 
radius  of  curvature  increases  to  a  value  R.  The  stresses  are  given  by 

<^0  «  Ey/R  (2) 

for  -n  <  y  <  n  and 

Oe  *  E[y/R  -  (yTn)/Ro]  (3) 

where  the  T  sign  pertains  to  the  regions  y  >  n  and  y  <  -p,  respectively. 

The  bending  moment  due  to  stress  distributions  of  Eqs.  (2)  and  (3),  which 
must  be  balanced  by  the  applied  moment,  is  given  by 

Ebh*  11  n  n 

M(R)  =  /  ogby  dy  = - [ - (1-3  -  +  4(-)>)]  (4) 

-h/2  12  R  Rq  h  h 

where  b  is  the  v<idth  of  the  beam.  Equations  (2)  through  (4)  can  also  be  used 
for  the  loading  cycle  by  setting  R  =  Rq  and  for  the  purely  elastic  case  by  addi¬ 
tionally  setting  n  =  h/2. 
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Equation  (4)  predicts  a  linear  dependence  of  1/R  on  the  applied  load  during 
unloading.  This  result  is  predicated  on  the  assumption  that  the  stress-strain 
relations  of  Eqs.  (2)  and  (3)  are  linear  (E  *  constant).  Thus,  experimental 
verification  of  Eq.  (4)  can  serve  as  a  check  to  which  extent  this  assumption  is 
warranted. 

In  analyzing  ultrasonic  velocity  changes,  knowledge  of  all  three  principal 
stresses  is  required.  However,  stresses  in  the  direction  perpendicular  to  the 
uniaxial  stress  can  be  neglected  here  since  they  are  expected  to  be  small. 
Boundary  conditions  require  these  stresses  to  be  zero  on  the  respective  perpen¬ 
dicular  surfaces.  Furthermore,  the  condition  of  plane  stress,  while  not 
strictly  applicable,  leads  to  the  compatibility  equation 


dCp  dOp 

«  Cp  +  r  -—  »  (jp  +  (R+y) 

dr  dy 


(5) 


Since  Cp  disappears  for  y  *  ±  h/2,  it  follows  from  Eq.  (5)  that  \  Op/OQ  1  ~  h/2R 
»  10"*  for  the  conditions  of  the  present  experiment. 

Velocity  Changes  Due  to  Elastic  and  Plastic  Deformation 

Hughes  and  Kelly  (ref  1)  have  derived  expressions  for  the  velocities  of 
ultrasonic  waves  in  elastically  stressed  solids  using  Murnaghan's  theory  of 
finite  deformations.  The  velocities  depend  on  the  triaxial  finite  strains 
through  the  Lame'  or  second-order  elastic  constants  X  and  u,  and  the  third-order 
constants  1,  m,  and  n. 

For  propagation  along  the  1-axis,  the  three  solutions  corresponding  to 
longitudinal  waves  and  transverse  waves  polarized  along  the  2-  and  3-axis, 
respectively,  are 


^0.  S.  Hughes  and  J.  L.  Kelly,  "Second  Order  Elastic  Deformation  of  Solids," 
Phvs.  Rev..  Vol.  92,  1953,  p,  1145. 
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PqVII  “  X  +  24  ♦  (21+X)0  +  (4m+4X+104)€i  (6) 

2  1 

PoVi2  *  M  +  (X+in)0  +  4»1£i  ♦  24€2  -  -  nes  (7) 

2  1 

Po'^ia  *  P  +  (X+m)0  ♦  44ei  +  2463  -  -  ne2  (8) 

Here,  Pg  is  the  initial  density  and  0  =  ej  +  £2  ■*■  £3- 

Johnson  (refs  2,3)  has  generalized  these  results  to  include  the  effects  of 
plastic  deformation  using  two  approaches.  His  first  approach  (ref  2)  leads  to 
the  addition  of  a  term  which  depends  on  the  plastic  deformations,  £^,  to  each  of 
the  Eqs.  (6)  through  (8).  The  resulting  velocities  may  be  written  as 


*  p 

PoVll.p  =  (X+24)4£i 

z 

*■  PoVll 

(9) 

PoVl2,p  *  4(2£i+2£2) 

+  PoVl2 

(10) 

PoVl3,p  =  4(2£i+2£§) 

+  PoVl3 

(11) 

Johnson's  second  approach  (ref  3)  leads  to  a  set  of  equations  which  are 
formally  identical  to  Eqs.  (6)  through  (8),  but  with  the  additional  proviso  that 
all  coefficients  in  these  equations  are  now  functions  of  the  plastic  defor¬ 
mation.  This  approach  then  predicts  that  velocity  changes  observed  during 
plastic  deformation  are  correlated  with  changes  in  the  second-  and  third-order 
elastic  constants.  In  both  approaches,  the  £.j  still  represent  the  elastic 
strains. 

If  we  make  the  additional  assumption  of  plastic  incompressibility,  Eqs.  (6) 
through  (11)  can  be  rewritten  in  terms  of  relative  velocity  changes. 

^G.  C.  Johnson,  "Acoustoelastic  Theory  for  Elastic-Plastic  Materials,"  J. 

Acoust.  Soc.  Am..  Vol .  70,  1981,  p.  591. 

^G.  C.  Johnson,  "The  Effect  of  Plastic  Deformation  on  the  Acoustoelastic 
Response  of  Metals,"  J.  Appl .  Mech..  Vol.  50,  1983,  p.  689. 


Specializing  for  the  case  of  uniaxial  strain  along  the  2-axis  and  setting  £^^3  = 
-y€2  «  -[X/(2X+2m) ]€2f  we  can  combine  the  various  theories  into  the  following 
set  of  equations 

^''11  u  mX  X* 

- »  oji  + . . (1-2X - )e2  (12) 

VqI  {X+2m){X+m)  u  ti 

^''12  1  nX 

- »  ai2  ♦ - (X+2u+m+  —  )C2  (13) 

Vqs  2(X+p)  4m 

^''13  1  X+M 

- «  0,3  ♦ . (m-2X - n)£2  (14) 

Vos  2(X>m)  2m 


Here,  Vqi  and  Vqs  are  the  longitudinal  and  shear  velocity  in  the  absence  of 
finite  stresses.  The  coefficients  o^j  are  zero  during  elastic  loading  and 
constant  for  elastic  unloading  after  plastic  deformation.  In  Johnson's  first 
approach  (ref  2),  we  have  ajj  =  013  «  -£2,  012  *  62/2.  and  all  the  elastic 
constants  remain  unchanged.  In  Johnson's  second  approach  (ref  3),  2ajj  = 
(AX+2Am)/(X+2m)  and  2oj2  »  20^3  »  Am/M-  1.  m,  and  n  can  also  change  during 
plastic  deformation. 

By  following  the  velocity  changes  through  the  elastic-plastic  and  then  the 
plastic-elastic  transition,  the  elastic  constants  1,  m,  and  n  can  be  determined 
before  and  after  plastic  deformation,  and  their  changes  can  be  correlated  via 
the  elastic  strain  £2  to  observed  velocity  changes  during  plastic  deformation. 
Thus,  some  definite  statements  about  the  various  theories  may  be  possible. 


2g.  C.  Johnson,  "Acoustoelastic  Theory  for  Elastic-Plastic  Materials," 
Acoust.  Soc.  Am.,  Vol .  70,  1981,  p.  591. 

^G.  C.  Johnson,  "The  Effect  of  Plastic  Deformation  on  the  Acoustoelastic 
Response  of  Metals,"  J.  Appl .  Mech.,  Vol.  50,  1983,  p.  689. 


EXPERIMENTAL  DETAILS 

The  specimen  used  was  a  36-inch  long  rectangular  bar  of  width  b  =  0.9077 
inch  and  height  h  «  1.95  inches,  made  from  ASTM  4340  steel  which  had  been  hot- 
rolled.  The  bar  was  austenitized,  quenched,  and  tempered  to  a  hardness  of  Rc 
31.  The  bar  was  subjected  to  four-point  bending.  The  support  and  the  load  were 
supplied  vertically  through  roller  bearings,  30-1/4  and  9-9/16  inches  apart, 
respectively.  Longitudinal  strain  was  measured  at  the  top  and  bottom  face  of 
the  bar  and  also  at  five  positions  evenly  spaced  across  the  height  of  the  bar. 
Transverse  strain  was  measured  at  the  top  and  bottom  face  of  the  bar.  Four 
1/4-inch  diameter  shear  transducers  were  mounted  with  their  centers  3/16-inch 
from  the  upper  and  lower  edge  of  the  bar,  with  polarization  along  the  length  and 
the  height  of  the  bar.  In  addition,  two  3/8-inch  diameter  longitudinal  trans¬ 
ducers  were  placed  with  their  centers  1/4-inch  from  the  top  and  bottom  edge. 
However,  only  one  of  the  longitudinal  transducers  retained  its  bond  throughout 
the  entire  experiment.  All  transducers  operated  at  5  MHz. 

The  load  was  typically  varied  in  1250-lb.  steps,  up  to  a  maximum  value  of 
21,250  lbs.  At  each  load,  strain  gage  readings  were  taken  and  corrected  for 
transverse  sensitivity.  Ultrasonic  time  increments  relative  to  the  unstressed 
state  were  measured  by  observing  the  time  shift  of  the  third  transverse  and  the 
sixth  longitudinal  echo  on  an  oscilloscope  using  the  20  nsec/cm  sweeprate. 

RESULTS  AND  DISCUSSION 

Figure  1  shows  the  load  versus  the  inverse  of  the  radius  of  curvature  as 
determined  from  a  linear  fit  to  the  longitudinal  strain  gage  readings.  The 
modulus  of  elasticity  derived  from  the  linear  region  during  loading  up  to  10,000 
lbs.  is  E  =  29.99  x  10*  lb. /in.*  ■  20.68  x  10*  bar.  This  is  in  close  agreement 
with  the  value  E  =  p(3X+2p)/(X+p)  =  20.61  x  10*  bar  obtained  using  X  =  11.03  x 
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10*  bar  and  m  »  7.99  x  10*  bar  from  an  earlier  experiment  on  a  similar  type 
steel  (ref  4).  Thus,  we  have  used  these  values  for  X  and  p  in  the  analysis 
below. 

As  can  also  be  seen  in  Figure  1,  the  transition  from  elastic  to  plastic 
deformation  starts  at  about  15,000  lbs.  During  unloading,  the  bar  promptly 
reverts  to  the  elastic  state  with  approximately  the  same  modulus  of  elasticity. 
However,  deviations  from  linearity  are  apparent  at  small  loads,  indicating  that 
there  is  reverse  yielding  of  favorably-oriented  grains.  Consequently,  we  have 
restricted  the  data  analysis  at  present  to  the  upper  part  of  the  unloading 
curve. 

The  load  at  minimum  radius  of  curvature  can  be  used  to  determine  the  loca¬ 
tion  of  the  elastic-plastic  interface  from  Eq.  (4).  We  find  |  n  1  “  0.40  in  good 
agreement  with  residual  stress  measurements  across  the  bent  bar. 

Relative  velocity  changes  were  determined  in  the  usual  manner  from  the 
relative  time  increments  and  the  transverse  strain,  calculated  at  the  detector 
position  from  the  strain  gage  readings.  They  are  shown  in  Figure  2  as  a  func¬ 
tion  of  longitudinal  strain.  The  result  for  the  shear  wave  polarized  along  the 
principal  stress  direction  (bottom)  is  particularly  noteworthy.  It  shows  a 
nearly  linear  dependence  on  strain  during  elastic  loading.  The  corresponding 
stress  scale  on  the  right-hand  side  is  obtained  using  the  value  of  E  given 
above.  At  17,500  lbs.,  a  sharp,  unusually  well-defined  change  in  slope  to 
nearly  constant  velocity  occurs,  indicating  a  drastic  change  in  material  behav¬ 
ior,  presumably  the  elastic-plastic  transition.  During  partial  unloading,  the 
material  reverts  promptly  to  the  elastic  state. 

*wT  Scholz  and  J.  Frankel,  "Acoustoelastic  Effects  in  Autof rettaged  Steel,"  in: 
Proceedinqs  of  the  1985  Ultrasonics  International  Conference.  Butterworth  & 
Co.,  Guilford,  Surrey,  UK,  1985,  pp.  441-446. 
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From  the  data  given  in  Figure  2  we  have  extracted  acoustoelastic  constants 
and  the  third-order  elastic  constants  1,  m,  and  n  with  the  help  of  Eqs.  (12) 
through  (14),  using  values  for  X  and  p  given  above  (ref  4).  In  deriving  these 
constants,  we  have  restricted  ourselves  to  the  central  region  of  the  elastic 
loading  and  the  top  part  of  the  elastic  unloading  curve.  Results  are  collected 
in  Table  I.  Since  the  systematic  errors  tend  to  be  self -compensating  in  the 
present  experiment,  the  differences  in  the  values  for  compression  and  tension 
and  before  and  after  plastic  deformation  are  believed  to  be  significant.  Stress 
acoustic  constants  can  be  obtained  from  these  acoustoelastic  constants  by 
dividing  by  E. 

TABLE  I.  ACOUSTOELASTIC  CONSTANTS  d(AVi j/Vq! ^s)/de2  AND  THIRD-ORDER  ELASTIC 
CONSTANTS  IN  10*  BAR. 


Compression 


Tension 


Unload 


Unload 


[2 

-1.641 

-1.444 

-1.753 

-1.815 

13 

0.128 

0.238 

0.000 

-0.102 

12-13 

-1.769 

-1.682 

-1.753 

-1.713 

LI 

0.210 

0.305 

0.140 

n 

-75.8 

-73.6 

-75.4 

-74.4 

m 

-63.3 

-56.5 

-67.7 

-70.4 

1 

-36.6 

-36.5 

-50.9 
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Table  I  depicts  the  acoustoelastic  constants  and  1,  m,  and  n  in  tension  and 


compression  as  well  as  changes  in  these  constants  with  plastic  deformation.  The 


least  variability  is  observed  in  n  which  is  the  only  third-order  constant  that 


enters  into  the  often  used  birefrinqent  expression.  Eg.  (13)  minus  Eg.  (141 


The  differences  in  the  corresponding  acoustoelastic  constant  (Table  I)  are  also 


rather  minimal. 


In  the  region  of  plastic  deformation,  all  velocities  appear  to  decrease 


with  the  amount  of  plastic  deformation.  From  the  intercepts  of  straight  line 


sections  fitted  to  the  data  points  in  Figure  2,  relative  velocity  changes  in  the 


stressed  state  can  be  estimated.  Using  Eq.  (3)  in  conjunction  with  the  value 


I  n  I  “  0.40  inch  for  the  elastic-plastic  interface,  the  strain  for  zero  stress  at 


the  transducer  position  can  also  be  estimated,  and  the  corresponding  relative 


velocity  changes  read  off  from  Figure  2.  Results  are  given  in  Table  II,  The 


observations  are  in  disagreement  in  magnitude,  as  well  as  in  relative  signs, 


with  the  predictions  of  Johnson's  first  approach  (ref  2)  as  discussed  in  con¬ 


junction  with  Eqs.  (12)  through  (14).  Furthermore,  they  cannot  be  explained 


solely  on  the  basis  of  changes  in  texture  since  for  any  texture, 


Po(vii*+vi2*+vi32 )  has  to  remain  constant,  and  decreases  in  some  velocities  must 


be  compensated  for  by  increases  in  others  (ref  5). 


In  Johnson's  second  approach  (ref  3),  relative  velocity  changes  during 


plastic  deformation  at  zero  stress  are  due  to  changes  in  X  and  pi,  and  the  dif¬ 


ferences  between  the  stressed  and  the  zero  stress  values  are  caused  by  changes 


‘G.  C.  Johnson,  "Acoustoelastic  Theory  for  Elastic-Plastic  Materials,"  J. 
Acoust.  Soc.  Am.,  Vol.  70,  1981,  p.  591. 

^G.  C.  Johnson,  "The  Effect  of  Plastic  Deformation  on  the  Acoustoelastic 
Response  of  Metals,"  J.  Appl.  Mech..  Vol.  50,  1983,  p.  689. 

^C.  M.  Sayers,  "Ultrasonic  Velocities  in  Anisotropic  Polycrystalline 
Aggregates,"  J.  Phys.  D;  Appl.  Phys..  Vol.  15,  1982,  p.  2157. 
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in  acoustoelastic  response,  i.e.,  changes  in  1,  m,  and  n.  Table  II  also  lists 
values  for  these  differences  calculated  from  Eqs.  (12)  through  (14)  with  the 
acoustoelastic  constants  from  Table  I. 


TABLE  II.  RELATIVE  VELOCITY  CHANGES  (AVij/Vos  l)  x  10»  DUE  TO  PLASTIC 
DEFORMATION.  EXPERIMENT,  UNDER  STRESS  (1);  EXPERIMENT,  ZERO 
STRESS  (2);  DIFFERENCE  (3)  «  (1)  -  (2);  CALCULATED  DIFFERENCE 
(4).  SEE  TEXT  FOR  DETAILS. 


i.i 

Compression 

Tension 

■n 

■n 

mm 

■n 

■n 

■n 

mm 

Hl^l 

(12) 

-1.52 

-0.17 

-1.35 

-0.88 

-0.83 

-0.89 

0.06 

O 

1 

(13) 

-0.65 

-0.51 

-0.14 

-0.49 

-1.95 

-1.52 

-0.43 

-0.44 

(12)-(13) 

-0.87 

0.34 

-1.21 

-0.39 

1.12 

0.63 

0.49 

0.17 

(11) 

-1.23 

-0.98 

-0.25 

-0.63 

Although  the  agreement  with  the  experimental  differences  is  not  perfect, 
the  numbers  are  remarkably  close  suggesting  that  there  is  some  basis  for 
Johnson's  second  approach  (ref  3).  It  must  also  be  realized  that  errors  of  the 
experimental  differences  may  be  fairly  large  and  that  residual  texture- induced 
velocity  changes  may  also  be  present.  Further  work,  experimental  as  well  as 
theoretical,  would  clearly  be  desirable. 


CONCLUSION 

The  present  experiment  demonstrates  that  four-point  bending  of  a  rec¬ 
tangular  beam  can  be  used  to  determine  the  acoustoelastic  response  of  ASTM  4340 
steel  during  elastic  and  plastic  deformation  in  both  tension  and  compression. 


^G.  C.  Johnson,  "The  Effect  of  Plastic  Deformation  on  the  Acoustoelastic 
Response  of  Metis,"  J.  Appl.  Mech,,  Vol,  50,  1983,  p.  689. 
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Third-order  elastic  constants  1,  m,  and  n,  derived  from  these  measurements,  show 
differences  between  tension  and  compression  and  changes  with  plastic  defor¬ 
mation.  Least  affected  is  the  constant  n  which  determines  the  acoustoelastic 
constant  for  birefringence  measurements.  The  measured  velocities  decreased  with 
plastic  deformation.  Measured  relative  velocity  changes  due  to  plastic  defor¬ 
mation  have  been  compared  with  predictions  of  various  theories. 
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